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Using a standard k-co turbulence model, in this paper detailed numerical computations 
have been carried out in inert simulators to examine the geometric dependence of transient 
flow features of solid rocket motors. We observed that the damping of the temperature 
fluctuation is faster in solid rocket with convergent port than with divergent port geometry. 
We inferred that an early damping of the transient flow fluctuations using the port geometry 
is a meaningful objective for the suppression and control of the instability and/or 
pressure/thrust oscillations during the staring transient of solid rockets. 


I. Introduction 

T HE basic idea behind a solid rocket motor (SRM) is simple but its design is a complex technological problem 
requiring expertise in diverse sub disciplines to address all of the physics involved. 1 " 21 The design optimization 
of the high-performance rockets is more complex when the mission demands dual-thrust. Many dual-thrust motors 
(DTMs) known to have experienced abnormal high ignition peak pressure often in the order of five times 
the steady state value during the ignition transient. 6,7,20,21 The motivation for the present study emanates from the 
desire to explain the phenomena or mechanism(s) responsible for the high ignition peak pressure (pressure spike), 
pressure-rise rate, instabilities and pressure oscillations often observed during the static tests and the actual flights 
of certain class of high-performance SRMs with non-uniform ports. 4 ' 9 Various measures were taken in the India 
industry to eliminate the peak pressure, but none of the conventional remedies seemed to help for the DTMs. 
Nevertheless, through the empirical techniques increasing the port area of the motor has been proposed as one of the 
remedies for reducing the unusual pressure spike of the DTM. Although such a remedy could negate the 
unacceptable peak pressure; it has affected the high-performance of the motor due to the envelop restriction. Hence 
the elimination of the unusual pressure spike and the pressure-rise rate without sacrificing the basic grain 
configuration or the volume loading became a meaningful objective for further studies. It may be noted here that 
none of the SRMs with convergent port geometry experienced abnormal high pressure spike and pressure-rise rate 
during the ignition transient. It has been shown conclusively through the previous experimental and numerical 
studies on SRMs with sudden expansion/convergence of ports that under certain conditions, secondary ignition can 
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occur far downstream of such ports. 6 " 8 One 
secondary ignition would results in two 
additional flame fronts, one spreading 
forward and the other backward. The 
effective time required for the complete 
propellant surface area to be ignited comes 
down drastically giving rise to high 
pressurization rate (dP/dt) in the second 
phase of ignition transient. This will lead to 
hard start of a rocket motor. Having logical 
precise relevance to the matter at hand, the 
next step through this paper will be to 
examine the intrinsic flow physics 
actuating the flow separation and 
secondary ignition in high-performance 
SRMs with non uniform port. Luke, G.D., 
et al., 9 reported that the U.S. Space 
Shuttle’s redesigned solid rocket motor 
(RSRM) head-end pressure rise rate is 
almost twice as high as the TITAN motors. 
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Figure 1. Comparison of the head-end pressure rise rate of the 
U.S. Space Shuttle’s redesigned solid rocket motor and various 
Titan motors (Configurations are shown in inset. Source: Ref. 9). 


The results from the existing models seem to be inconclusive as to the cause of the relatively higher-pressure rise 
rate encountered in the Space Shuttle motor (RSRM). Note that the head-end port of the RSRM is narrow, compared 
to the TITAN motors. Moreover RSRM port geometry is having divergent region, whereas Titan motors are having 
port with convergent region (see Fig. 1). In this paper using the inert simulators of SRM, parametric studies have 
been carried out to examine the influence of port geometries on instabilities and pressure oscillations during the 
starting transient. 


II. Overview of the Numerical Methodology 

Numerical simulations have been carried out with the help of a well-established two-dimensional standard 
k-omega model. This turbulence model is an empirical model based on model transport equations for the turbulence 
kinetic energy and a specific dissipation rate. This code solves standard k-omega turbulence equations with shear 
flow corrections using a coupled second order implicit unsteady formulation. This model uses a control-volume 
based technique to convert the governing equations to algebraic equations. The viscosity is determined from the 
Sutherland formula. An algebraic grid-generation technique is employed to discretize the computational domain. A 
typical grid system in the computational region is selected after the detailed grid refinement exercises. The grids are 
clustered near the solid walls using suitable stretching functions. The motor geometric variables and material 
properties are known a priori. Initial wall temperature, inlet total pressure and temperature are specified. At the solid 
walls a no-slip boundary condition is imposed. At the nozzle exit a pressure profile is imposed. The Courant- 
Friedrichs-Lewy number is initially chosen as 3.0 in all of the computations. Ideal gas is selected as the working 
fluid. The transient mass addition due to propellant burning is deliberately suppressed in this model to examine the 
turbulent separated flow features discretely in SRMs with divergent ports. The code has successfully validated with 
the help of benchmark solutions. 7,8,18 


III. Results and Discussion 

In the first phase of this study low-velocity transient (LVT) motors (A t /A p < 0.56, L/D < 10) and in the second 
phase high-velocity transient (HVT) motors (A t /A p > 0.56, L/D >10) with non-uniform ports are considered. In 
order to avoid the complications arising from propellant combustion, in both cases, inert simulators are considered 
for examining the internal flow features of SRMs with non-uniform ports. Figure 2 is demonstrating the influence of 
the port geometry on the axial velocity variations of five different LVT motors with different port geometry but with 
same initial and boundary conditions. In the first three cases divergent location (X s ) is varied, in the forth case inlet 
diameter is increased by 50% and in the fifth case divergence angle, a is increased from 45° to 64°. All the results 
reported are anticipated and giving corroborative evidences of the previous experimental and theoretical findings. 6 ' 8 
It can be seen from Fig.2 that, in three different cases (Case 2, 3 & 5), the axial velocity is relatively high 
near the divergence location. This can be explained with the help of boundary layer theory. Note that owing to the 
viscous friction, boundary layer will be formed on the walls (before the transition region) and their thickness will 
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increase in the downstream direction to the 
divergence location. Since the volume of flow 
must be the same for every section, the decrease 
in rate of flow near the walls which is due to 
friction must be compensated by a corresponding 
increase near the axis. Thus the boundary layer 
growth occurs under the influence of an 
accelerated external flow. As a result, at larger 
distances from the inlet section velocity will be 
relatively high and the flow will possibly become 
turbulent; consequently the boundary layer 
thickness will suddenly increase leading to the 
sudden increase in the axial velocity due to the 
rocket motor port area fraction blocked by the 
boundary layer displacement thickness. This will 
cause flow separation far downstream of the 
divergence region. In the fourth case reported 
herein shown relatively low velocity at the axis 
due to the high port area compared to the other 
four cases reported. As stated in the introduction 
increasing the port area of a solid rocket motor 



Figure 2. Velocity variations along the axis of LVT motors 
showing the peak near the transition location (Baseline 
values: L/D = 4, X s /d = 3, d t /D = 0.375). 


can reduce the unacceptable peak pressure, caused due to the gas dynamics of the upstream narrow port, at the 
expense of propellant loading density. 

In the first case (Xs = 0.03m, see Fig. 2) flow recirculation tendency, leading to reattachment and secondary 
ignition, was found very less because location of the transition region was near to the head-end at the forfeit of the 
propellant loading density. When the transition location was fixed at far downstream of the SRM, the tendency of 
flow separation was found very high. This will lead to the formation of recirculation bubble and flow reattachment. 
Note that the flow reattachment will favor secondary ignition and that will cause the flow unsteadiness leading to an 
unacceptable high-pressure rise rate during the staring transient period of operation of solid rockets. Hence the 
prudent selection of the transition location within the given envelop, without diluting the high-performance nature of 
solid rocket motor, is critical for a designer. This task will be more complex in the case of HVT motor, which is 
discussed in the subsequent session. 

In the second phase; SRMs with three different transition locations (X s ) are considered. Figure 3 is 
demonstrating the difference in velocity magnitude along the axis of HVT motors with different port geometries at 

two different time intervals but with same initial 
and boundary conditions. In all the cases the 
velocity magnitude is found maximum at the 
transition location. As explained in the previous 
cases, in general, at larger distances from the 
inlet section (X s ) velocity will be high at the 
step location due to the boundary layer effect. 
But Fig.3 shows that when Xs = 0.20 m the 
peak value of the axial velocity is relatively 
lower than the former case (Xs = 0.15 m). This 
will not contradict the argument reported earlier 
based on boundary layer effect because this 
difference is due to the altered variation of 
the entire flow field due to the nozzle end 
effect coupled with the geometry dependent 
driving forces and the corresponding 
compressibility effect. Through these diagnostic 
investigations, we observed that there is a 
limiting distance of X s for developing 
maximum axial velocity in any HVT motor due 
to its port geometry coupled with nozzle end 
effect. In all the HVT motor cases, as 
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Figure 3. Velocity variations along the axis of three HVT 
motors showing that the high peak is near the transition region 
(A typical grid system is shown in inset). 
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anticipated, at the upstream the flow acceleration is found very high compared to the LVT motor cases reported 
earlier. 

Note that the thickness of a turbulent boundary layer is larger than that of a laminar boundary layer owing to 
greater energy loss in the former. The development of the wall boundary layer in turbulent flow is more complicated 
than in wholly laminar flow. Initially it takes the form of a laminar layer, but at some position along the rocket 
motor port there is a transition to a turbulent layer, where a sudden increase in axial velocity can be discerned (see 
Fig.2, Cases 2, 3 & 5). The actual position of transition depends on a number of factors including Reynolds number, 
surface roughness, and the turbulence level of the igniter jet flow entering the motor port. 
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Figure 4. Shifting of the reattachment point towards 
the step at high-turbulence level (Enlarged view at the 
transition region of an HVT motor). 
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Figure 5. Demonstrating the variation of the 
turbulence intensity along the axial direction of 
an HVT motor; and pinpointing that the 
turbulence intensity is maximum near the eye of 
the recirculation bubble (Velocity vectors are 
shown in inset). 


In another attempt the influence of the igniter jet turbulence intensity (JTI) on flow separation has been 
examined. A case with high-turbulence intensity shows less possibility of flow separation and reattachment. Figure 4 
shows the near wall velocity viewing the tendency of the shifting of the reattachment point towards the step location 
at higher turbulence level. We also observed that in all the cases the maximum turbulence intensity falls at the eye 
(not the center point!) of the recirculation bubble (see Fig. 5) because within the separation bubble, the mean 
turbulent intensity rises highly especially near the centre of the bubble due to high mixing and large-scale 
unsteadiness, and it reduces as one move to the reattachment point and over the boundary layer development region. 

Note that the secondary ignition occurs inside the initial recirculation bubble. Therefore, the size of the 
recirculation bubble and the location of the eye are important for further study. It is inferred that by inducing high 
turbulence level (practically by using igniter nozzle wire screens, or by creating artificial roughness to the grain wall 
without affecting the ballistics!), the position of transition could be brought closer to the entry region, or indeed the 
laminar layer could be entirely eliminated! Separation is mostly an undesirable phenomenon because it entails large 
energy losses. 



Figure 6. A typical HVT motor with divergent port. 



Figure 7. A typical HVT motor with convergent port. 
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In another attempt flow features of an HVT motor with divergent port geometry has been compared with 
convergent port under identical inflow conditions. The configurations of these motors and the corresponding grid 
systems in the computational domain are shown in Fig. 6 and Fig. 7 respectively. We observed that at the identical 
inflow conditions the magnitude of the velocity fluctuation is more dominant for the rocket motor with divergent 
port than with the convergent port (see Figs. 8 & 9). It is more prominent at the transition region of the grain (see 
Fig. 8). Figure 10 shows the comparison of the dynamic pressure along the axis of convergent and divergent port 
cases. It is evident from this figure that at the same inflow conditions the magnitude of the dynamic pressure is 
relatively high in the case of an HVT motor with divergent port. This is more pronounced at the transition region. 



Axial Distance, m 

Figure 8. Demonstrating the axial velocity 
variations of an HVT motor with divergent port at 
different time intervals (configuration is shown 
along with the x-axis). 
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Figure 9. Demonstrating the axial velocity 
variations of an HVT motor with convergent port 
at different time intervals (configuration is shown 
along with the x-axis). 



Figure 10. Comparison of the dynamic pressure 
variations at the axis of two HVT motors at different 
time intervals and with same initial and boundary 
conditions. 



Figure 11. Comparison of the total temperature 
variations at the axis of two HVT motors at 
different time intervals but with same initial and 
boundary conditions. 


Figure 11 shows the comparison of the total temperature at the axis of two different HVT motors at four different 
time intervals. It is evident from this figure that the steady state temperature will attain rapidly in the case of HVT 
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motor with convergent port than with divergent 
port. It can also be seen that the temperature 
fluctuation is more prominent for divergent 
port case due to the flow unsteadiness 
associated with separation and reattachment. 
Figures 12-15 show the contours of static 
pressure, dynamic pressure and total 
temperature at two different time intervals for 
the divergent and convergent port cases 
respectively. When time advances the flow 
field significantly altered in both cases. The 
recirculation bubble is discemable in the case 
of divergent port case at time t = 0.0005 s and it 
further expands downstream when time 
advances. These changes are reflected in the 
dynamic pressure contours of the HVT motor 
with divergent port (see Figs. 12 &13). Note 
that the downstream of the separation point, the 
velocity u close to the surface becomes 
negative and so a region of reverse flow is 
established. Because of their ability to transfer 
momentum laterally, turbulent flows are more 
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Figure 12. Demonstrating the critical flow features of an 
HVT motor with divergent port at an early period of starting 
transient (t = 0.001 s). 

able than laminar flows to negotiate regions of adverse pressure gradients. Whether or not separation actually takes 
place, the general effect of the adverse pressure gradient is to give rise to a localized region of slow moving fluid 
stretching away from the wall. Because of the continuity condition, which can be applied over the whole cross- 
sectional area, the axial flow velocities must necessarily increase elsewhere in order to compensate for this effect. 
There is therefore a tendency for flows to become increasingly non-uniform whenever positive axial pressure 
gradients are encountered. 

From these flow features (see Figs.8-15), the rocket motor designer can ascertain that there is a possibility of 
more thrust oscillations in HVT motors with divergent port than with convergent port geometry. Although such 
oscillations do not peril the launcher operation, they induce some penalties to the overall performance. Hence, in the 
mission point of view, to avoid the unpredictable performance variation one has to recast the intrinsic shape of the 
grain geometry at the transition region in correspondence with the mechanical response of the propellant based on 

the magnitude of the anticipated dynamic 
pressure/thermal loads acting on the 
propellant surface. 1 The suppression or the 
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1.49e+05l meaningful objective for any solid rocket 
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Figure 13. Demonstrating the critical flow features of an 
HVT motor with divergent port at a later period of starting 
transient (t = 0.004 s). 


motor designer. While comparing the 
characteristic curves, one can easily notice 
that under the identical inflow conditions a 
case with divergent port shows high 
magnitude of velocity along the axis. 
Unlike the divergent case, the total 
temperature variation for the convergent 
port case is found amplifying along the axis 
during the initial period like a sinusoidal 
wave (see Fig. 11). Further the magnitude 
of the total temperature got diminished and 
reached the steady state value faster than 
the divergent port case presumably for 
readjusting it with the corresponding low 
velocity and pressure values, which has 
altered due to the port geometry. It was 
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seen that the igniter temperature fluctuations 
are diminished rapidly in the case of 
convergent port case irrespective of the 
igniter total pressure. Note that the total 
temperature has to be readjusted with the 
corresponding velocity and the pressure 
values within the motor, which has altered 
due to the sudden variation of the port 
geometry. From these findings we can 
anticipate a complex physical situation in 
HVT motors with non-uniform port during 
the transient combustion, which however is 
not considered in this analysis. These 
parametric studies lead to say that a prudent 
selection of the port geometry is one of the 
challenging tasks for any HVT motor Figure 14. Demonstrating the critical flow features of an 
designer. HVT motor with convergent port at an early period of starting 

Note that near to a solid surface flow transient (t = 0.001 s). 
velocities are low due to the no-slip condition 

at the wall. Hence in a region where the piezometric pressure is increasing, there are likely to exist certain 
streamlines, on which there are points whose total pressure is less than the piezometric pressure a little further 
downstream. When this happens, these streamlines can only reach this further point if their energy is increased by 
the action of the shear force exerted by adjacent elements of the flow. This condition is satisfied when dx/d y>0, 
where x is the local shear stress and y is the distance measured away from the grain wall. This process of energy 
conversion by the action of viscosity cannot be maintained indefinitely and, if the flow does not manage to negotiate 
the region of adverse pressure gradient, a point is reached at which the value of x and hence of du/dy becomes zero 
at the surface. Downstream of such a point, which is known as a separation point, the velocity u close to the surface 
becomes negative and so a region of reverse flow is established. Because of their ability to transfer momentum 
laterally, turbulent flows are more able than laminar flows to negotiate regions of adverse pressure gradients. 
Whether or not separation actually takes place, the general effect of the adverse pressure gradient is to give rise to a 
localized region of slow moving fluid stretching away from the wall. Because of the continuity condition, which can 
be applied over the whole cross-sectional area, the axial flow velocities must necessarily increase elsewhere in order 
to compensate for this effect. There is therefore a tendency for flows to become increasingly non-uniform whenever 
positive axial pressure gradients are encountered (see Figs. 12-15). 

We also discerned that at the subsonic inflow conditions, the flow gets accelerated to a higher Mach number 
(M>1) near the transition region of a typical HVT motor having narrow upstream port ending with divergent port 

but without any geometrical-throat! A shock 
wave cannot exist unless the Mach number is 
supersonic; therefore the flow must have 
accelerated through a throat which is sonic. 
This is presumably due to the formation of 
the fluid-throat at the beginning of the 
transition region due to the area fraction 
blocked by the boundary layer displacement 
thickness. As a result the upstream narrow 
port of the DTM could act like a second 
igniter to the downstream port leading to the 
formation of the possible shock waves inside 
the motor during the starting transient period. 
Note that the downstream of the shock the 
Figure 15. Demonstrating the critical flow features of an flow has an adverse pressure gradient, usually 
HVT motor with convergent port at a later period of starting leading to wall boundary-layer separation and 
transient (t = 0.004 s). reattachment. The shock waves in HVT 

motor will alter the turbulence level and this 
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new turbulence level will alter the location of reattachment and secondary ignition. These are discussed in the 
companion paper. 21 

From these studies one can deduce that the thrust/pressure oscillations, pressure-rise rate and unexpected peak 
pressure often observed in SRMs with non-uniform ports are presumably contributed due the joint effects of the 
geometry dependent driving forces and the chamber gas dynamic forces. The present study is expected to aid the 
designer for conceiving the physical insight into problems associated with the prediction and the reduction of the 
peak pressure, the pressurization rate and thrust oscillations during the starting transient period of operation of 
SRMs with non-uniform ports. 


IV. Concluding Remarks 

The fact that a separated flow region is formed downstream of a sudden expansion area is easy to appreciate. 
Rather less obvious is the fact the flow may separate from a surface which has no discontinuities of curvature. The 
process of separation is associated with large rates of dissipation of mechanical energy and so the avoidance of 
separation is an important factor in the design of solid rockets. We concluded that due to the area fraction blocked 
by the boundary layer displacement thickness the SRMs axial velocity will be increasing up to a critical distance 
from the inlet section and further it will alter depends upon the port geometry and the nozzle-end effect, which in 
turn alter the separation and secondary ignition. We inferred that the location of the reattachment and/or secondary 
ignition point and the features of the recirculation bubbles are sensitive to the igniter jet turbulence intensity. We 
also inferred that under identical inflow conditions head-end chamber pressurization-rate (dP/dt) is found relatively 
less in the case of SRMs with convergent port geometry than with divergent cases. We discerned that at the subsonic 
inflow conditions there is a possibility of the occurrence of shock waves inside the port of a rocket motor due to the 
formation of the fluid-throat as a result of the area fraction blocked by the boundary layer displacement thickness. 
The shock waves will further alter the location of the reattachment/secondary ignition point and also enhance the 
heat flux to the propellant surface, which in turn enhance the flame spread rate and the transient burning. This will 
lead to the unexpected peak pressure, the pressure-rise rate and the thrust oscillations during the mission of HVT 
motors with non-uniform port. We also discerned that the situation will be more complex during the transient 
combustion, which however is not considered in this study. It was seen that the igniter temperature fluctuations are 
diminished rapidly in the case of convergent port case irrespective of the igniter total pressure. Note that the total 
temperature has to be readjusted with the corresponding velocity and the pressure values within the motor, which 
has altered due to the sudden variation of the port geometry. We concluded that the damping of the total temperature 
fluctuation within the motor is a meaningful objective for the suppression and control of instability and/or 
pressure/thrust oscillations often observed in HVT motors with non-uniform port. These parametric studies lead to 
say that a prudent selection of the port geometry is one of the challenging tasks for any dual-thrust motor designer. 
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